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ABSTRACT 

High resolution spectra of the hot central star of the planetary nebula (CSPN) M27, acquired with the Far 
Ultraviolet Spectroscopic Explorer (FUSE), have revealed an unusually rich set of narrow molecular hydrogen 
absorption features. This object is also unique in that the velocity of nebular absorption features are cleanly 
separated from the velocity of the intervening interstellar medium. These features blend with and in many cases 
obscure atomic features. We have developed a continuum model of the CSPN including atomic and molecular 
hydrogen absorption. Using this model we have identified and tabulated the metal lines as arising from either 
photospheric, nebular and/or non-nebular velocity systems. We find the nebular outflow and ionization balance 
to be stratified with high ionization states favored at low velocity and low ionization states favored at high 
velocity. Neutrals and molecules are found at a velocity that marks the transistion between these two regimes. 
These observations are a challenge to the interacting wind model of PN evolution. Mappings at high resolution 
of the line profiles for C I - IV, N I - III, O I, O VI, Si II - IV, P II - V, S II - IV, Ar I - II and Fe II - III within the 
FUSE and STIS bandpasses are presented. The digitial spectra of the star and the model are freely available 
on the H200IS website. They will be useful for photospheric analyses seeking to determine the metallicity of 
the central star and for absorption line based atomic and molecular abundance determination of the nebular 
outflow. 

Subject headings: atomic processes — ISM: abundances — (ISM:) dust, extinction — (ISM:) planetary nebu- 
lae: general — (ISM:) planetary nebulae: individual (NGC 6853) — line: identification — 
line: profiles — molecular processes — plasmas — (stars:) circumstellar matter — (stars:) 
white dwarfs — ultraviolet: ISM — ultraviolet: stars 



I. INTRODUCTION 

High resolution spectra of the hot central star (CS) of the 
Planetary Nebula (PN) M27 (NGC 6853, PN G060.8-03.6, the 
Dumbbell), acquired with the Far Ultraviolet Spectroscopic 
Explorer (FUSE), have revealed an unusually rich set of nar- 
row molecular hydrogen absorption features, spanning the en- 
tire bandpass. FUSE carriers no on-board source for wave- 
length calibrations and consequently M27 has been observed 
numerous times for this purpose, producing in the process a 
high quality data set as illustrated in Figure [T] The high den- 
sity of the molecular lines presents problems for the identifi- 
cation and analysis of atomic features arising from the stellar 
photosphere and inters tellar medium (ISM). 

In a separate work (iMcCandliss et al.ll2007L hereafter Pa- 
per I) we have developed a molecular and atomic hydro- 
gen absorption model that, when applied to a readily avail- 
able pure hy drogen helium hot white dwarf model of the CS 
dRauchl 1200 3). provides a reasonable template against which 
flie atomic absorption features from other metal line systems 
can be isolated. Here we compare the model to the obser- 
vations and present a list of line identifications for atomic 
species associated with the stellar photosphere, the surround- 
ing nebula and non-nebular ISM. We also present high reso- 
lution absorption line profiles from STIS and FUSE for the 
metals C I - IV, N I - III, O I, O VI, Si II - IV, P II - V, S II - IV, 
Ar I - II and Fe II - III. 

Paper I finds the observed atomic and molecular velocity 
stratification in the nebular outflow of M27 is challenging to 
explain in the context of the standard interacting winds model 
for PNe (iKwok et al.l 119781) . A similar challenge has been 
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made by iMeaburnI ( l2005h (see also iMeaburn et al.ll2005allbl) 
who used high resolution position-and-velocity spectroscopy 
of the optical emission lines to derive the ionization kinemat- 
ics of several objects. They find that ballistic ejection could 
have been more important than interacting winds in shap- 
ing the dynamics of PNe. In the interacting winds model a 
high speed (~ 1000 km s^^) highly ionized radiation driven 
wind from the hot star shocks and ionizes a slow moving (~ 
10 km s^'), high density mostly molecular asymptotic giant 
branch (AGB) wind. In Paper I no evidence was found for a 
high speed radiation driven wind. The upper limit to the ter- 
minal expansion velocity is « 70 km s^^ The atomic line 
profiles show that high ionization material appears at low ve- 
locities (< 33 km s^') and low ionization material appears 
at high velocities (33 - 65 km s^^). At the transition veloc- 
ity (33 km s^') between these two regimes, a predominately 
neutral atomic medium appears along with the vibrationally 
excited H2. The molecular abudance in the nebula is low 
(A?(H2)/A?(HI)< 1). 

The Paper I study found that, despite the close proximity of 
a hot central star, the molecular hydrogen ro-vibration levels 
of the ground state were close to thermal ('^ 2000 K). They 
found no detectable levels of far-UV fluorescence of molec- 
ular hydrog en, although Lyg pumped fluorescence has been 
detected by iLupu et al.l (l2006h . Surprisingly, they also found 
that the line-of-sight was devoid of dust. The temperature im- 
plied by the molecular hydrogen suggests that the diffuse neb- 
ular medium is too hot for dust to form. They argue that the 
clumpy medium seen in the CO map of Bachiller et al. (2000), 
which cast the [O III] shadow s seen by Meaburn & LoDe3 
( 119931) and lO'Dell etal.1 ( |2002|) . are the source of the excited 
molecular hydrogen for the diffuse nebular medium and may 
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TABLE 1 
FUSE Observation Summary 



Observation 


Date 


APER 


EXPiLiFla 


EXP:SiClb 


EXP:LiF2a 


EXP:SiC2a 








(ksec) 


(ksec) 


(ksec) 


(ksec) 


Ml 070101 


2000-09-01 


MDRS 


15.5 


8.0 


15.6 


8.9 


M 1 07010? 


2000-09-24 


HIRS 


17.5 


4.8 


13.3 


5.6 


M 1 070101 


2001 -05-2X 


LWRS 


7.4 


7.5 


7.3 


7.0 


M 1 070104 


2001 -05-29 


MDRS 


7.0 


4.8 


6.6 


0.6 


M 1 070105 


2001 -05-29 


HIRS 


8.0 


4.1 


6.3 


2.0 


ivi i\j / \jj\jyj 


9001 -07-9S 


LWRS 


7.1 


7.1 


7 9 


7.1 


M 1 070^07 

iVl 1\J 1 \JJ)\J 1 


9001 -07-9S 


MDRS 


6.7 


3.5 


6.5 


5.1 


iVl 1\J 1 \JJ)\JO 


9001 -07-9Q 


HIRS 


8.5 


3.4 


6.3 


3.3 


Ml 070109 


2001-08-01 


LWRS 


1.5 


7!6 


1.1 


7!4 


Ml 0701 10 


2001-08-01 


MDRS 


6.6 


4.1 


6.4 


4.8 


Ml 0701 11 


2001-08-01 


HIRS 


7.8 


0.6 


5.6 


0.8 


Ml 0701 12 


2002-10-10 


LWRS 


5.0 


5.1 


4.8 


5.0 


Ml 0701 11 


2002-10-10 


MDRS 


7.0 


3.3 


6.4 


3.0 


M1070114 


2002-11-02 


HIRS 


3.3 


0.9 


2.4 


0.7 


Ml 0701 15 


2004-05-25 


LWRS 


5.3 


5.3 


4.3 


4.4 


Ml 0701 16 


2004-05-25 


MDRS 


4.0 


2.9 


4.9 


4.4 


Ml 0701 17 


2004-05-25 


HIRS 


6.5 


3.2 


5.0 


4.3 


Ml 0701 19 


2002-11-01 


HIRS 


2.7 


1.1 


2.0 


0.9 


PI 041101 


2000-06-05 


LWRS 


16.9 


16.5 


16.3 


16.9 



harbor depleted material (dust). A depletion analysis of the 
Fe, Si, and S in the nebular outflow could provide information 
as to whether these clumps are composed of depleted material 
that is being released into the diffuse medium along with the 
molecular hydrogen. 

The digitial spectra used in this study, are available at the 
following url, 

http : //www.pha . jhu . edu/"'stephan/H2ools/M27kruk04/. 
This dataset will be useful for photospheric analyses of the 
central star and for studies seeking to understand the kine- 
matics of the atomic and molecular abundances in the nebular 
outflow. 

2. DATASETS 

2.1. FUSE 

FUSE spectra of the CS (MAST object ID GCRV12336) 
were acquired under the observatory wavelength calibration 
program (ID: M107). The CS was observed numerous times 
through all three focal-plane apertures: large (30")^, medium 
(20" X 4") and small (20" x 1.25"); referred to respectively 
as LWRS, MDRS, and HIRS. A summary of the observa- 
tions is given in Table [T] All spectra were acquired in TTAG 
mode. Each observation was performed as a series of expo- 
sures, ranging in length from 0.5 ksec to 3 ksec. Each expo- 
sure produces 8 individual spectra; one from each of two de- 
tector segments for each of the 4 optical channels. Effective 
exposure times are listed in the table for the more interesting 
segment in each channel. Exposure times occasionally dif- 
fer slightly from one segment to another because the detector 
high voltage in a given segment can be reduced temporarily 
in response to event bursts or brief current spikes. Exposure 
times also differ from one channel to another as a result of 
channel misalignment effects. For a description of the de- 
tectors, channel alignment issues, system s nomenclature and 
oth er aspects of the FUS E instrument, see iMoos et aH (l2000h 
and lSahnow et al.1 (l2000l) . 

The exposures were processed with the following proce- 
dure. The raw data were calibrated using Calfuse 3.0.7. 
The resulting spectra for a given aperture (LWRS, MDRS 
or HIRS) were screened to identify the peak detected stellar 
flux. Exposures with less than 40% of the peak flux (a small 



fraction) were discarded. The remaining were normalized to 
the peak flux, yielding an effective exposure time. Chan- 
nel/segments were coaligned with a cross-correlation proce- 
dure on narrow ISM lines. Spectra were shifted by whole 
pixel units to avoid smoothing and combined by effective 
exposure time weighted average. At this point there is one 
spectrum for each spectrograph aperture, divided among the 
8 channel/segment combinations. The MDRS spectra were 
put on an absolute wavelength scale by adjusting each chan- 
nel/segment spectrum so that the H2 lines arising from the hot 
component o f the molecular gas are at a fiducial offset (-69 
km s^\ see § 14.1 [ regarding this point). Finally, the normaliza- 
tion of the MDRS and HIRS were readjusted to match that of 
flie LWRS and the zero point offsets of the HIRS and LWRS 
channel/segment were adjusted to match the MDRS.' We use 
the HIRS data in this work as it yields the highest resolution. 

Each extracted segment has associated with it a set of one 
dimensional arrays: wavelength (A), flux (ergs cm^^ s^' 
A^'), and estimated statistical error (in flux units). Figure [T] 
shows the signal-to-noise and spectral range for each of the 
coadded HIRS segments, obtained by dividing a flux array by 
the corresponding statistical error array. The signal-to-noise 
estimated in this way is purely statistical and does not account 
for systemic errors, such as detector fixed pattern noise. In the 
low sensitivity SiC channels the continuum signal-to-noise 
ranges between 10 - 25, while in the higher sensitivity LiF 
channels it is 20-55. The resolution of the coadded spectra 
changes slightly as a function of wavelength for each segment. 
In modeling the molecular hydrogen absorption, we find that 
a gaussian convolution kernel with a full with at half max- 
imum of 0.05 A at 1000 A (spectral resolution R = 18,000, 
velocity resolution bV « 17 km s^') provides a good match 
to the unresolved absorption features throughout much of the 
bandpass. 

' After preforming these procedures spectral overlays the HIRS LiF lb and 
LiF2a segments were in obvious disagreement with the LWRS and MDRS 
LiFlb and LiF2a segments. Consequently the HIRS wavelength scale was 
been adjusted by adding in a linear dilation of the scale of 0.01 A for every 
45 A starting at 1 125 A. This adjustment to the HIRS wavelength scale has 
been incorporated into subsequent versions of CALFUSE. The HIRS LiF2a 
flux was also rescaled by a factor or 0.91 to achieve agreement with LiFlb. 
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Signal-to-noise ratio of tlie eight coadded spectral segments taken through the HIRS spectrograph aperture. The individual segments are labeled. 



Close comparison of the wavelength registration for over- 
lapping segments reveals isolated regions, a few A in length, 
of slight spectral mismatch (~ a fraction of a resolution ele- 
ment) in the wavelength solutions. Consequently, combining 
all the spectra into one master spectrum will result in a loss 
of resolution. However, treating each channel/segment indi- 
vidually increases the bookkeeping associated with the data 
analysis. Further, because fixed pattern noise tends to domi- 
nate when the signal-to-noise is high, there is little additional 
information to be gained in analyzing a low signal-to-noise 
data set when high signal-to-noise is available. For these rea- 
sons we elected to form two spectra each of which covers the 
900 - 1190 bandpass contiguously, using the following pro- 
cedure. The flux and error arrays for the LiFla and LiFlb 
segments were interpolated onto a common linear wavelength 
scale with a 0.013 A bin, covering 900 - 1 190 A. The empty 
wavelength regions, being most of the short wavelength re- 
gion from 1000 A down to the 900 A and the short gap re- 
gion in between LiFla and LiFlb, were filled in with most 
of SiC2a and a small portion of SiC2b respectively. We re- 
fer to this spectrum as sl2. The LiF2b, LiF2a, SiCla and 
SiClb segments were merged similarily into a spectrum, s21. 
Absorption line analyses were carried out using both sl2 and 
s21. 

2.2. STIS 

The E140M spectrum of M27 (o64d07020_xld.fits), taken 
from the Multimission Archive at Space Telescope (MAST), 
was acquired for HST Proposal 8638 (Klaus Werner - PI). 
These data are a high level data product consisting of arrays 
of flux and wavelength calibrated, one-dimensional extrac- 
tions of individual echelle orders. The data were acquired 



through the 0"2x0"2 spectroscopic aperture for an exposure 
time of 2906 s and were reduced with CALSTIS version 2. 18. 
The spectral resolution of t he E140M is g iven in the STIS 
data handbook (version 7.0 lOuii and 120031) as w 48,000 
(5y = 6.25 km s- 1 ). The intrinsic line profile for the 0"2 x 0;'2 
aperture has a gaussian core with this width, but there is non- 
negligible power in the wings. These data offer a high velocity 
resolution view of the ionization stratification in the nebular 
expansion. 

3. MODEL OF CONTINUUM ABSORPTION BY HYDROGEN AND 
LINE IDENTIFICATIONS 

A model of the CS SED, including the atomic and 
molecular hydrogen absorptions has been developed by 
McCandliss et al. (2007) to aid identification of photospheric, 
nebular or non-nebular velocity components in the FUSE and 
STIS metal line profiles. We used synt hetic stellar flux inter- 
polated from the grid of lRauchI (|2003|) with \ogg = 6.5, T = 
120,000 K, and ratio of H/He = 10/3 by mass. This model in- 
cludes no metals and is consistent with, although slightly hot- 
ter than, the quantitative spectroscopy of iNapiwotzkil (I1999I) . 
The adopted temperature and gravity has slightly less pres- 
sure broadening and is a better match to the observed Lyman 
lines towards the series limit. We also adopt a stellar mass of 
0.56 Mq as suggested by post asymptotic giant branch evolu- 
tionary tracks. Use of this mass along with the above gravity 
required a distance of 466 pc t o match the absolute fl ux, which 
is at the upper limit given by [Benedict et alj ( l2003h . Beyond 
our immediate need to match the absolute flux for the given 
gravity and mass there is no particular reason to prefer this 
distance over that derived by Benedict et al. Questions re- 
garding the acceptable uncertainty in distance, absolute flux 
and derived stellar parameters are best left for a stellar model 
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Fig. 2a. — Hydrogen absorbed continuum model 905 - 980 A. Identified line blends are numbered and marked with upper and lower + and tabulated in Table 
12] Upper portion of each panel is an overplot of the sl2 (orange) and s21 (green) FUSE spectra, and the model (black). Lower portion of each panel shows the 
data after subtraction of the model, offset by 10^'- ergs cm"- s"' A"'. 
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Fig. 2b. — Hydrogen absorbed continuum model 975 - 1050 A. See Figure 2a for a description of the colors and symbols. 
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Fig. 2c. — Hydrogen absorbed continuum model 1045 - 1 120 A. See Figure 2a for a description of the colors and symbols. 
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specifically tailored to include the effects of metals, gravity 
and evolutionary state. The data supplied here is intended to 
enable such an effort. 



TABLE 2 

Stellar, Nebular and Non-Nebular Line ID'S 





A 


Atnmir TO ^ 
A 







1171.85 


O VI 1172.439? 


Ph? 


1 


1170.95 


O VI 1171.561? 


Ph? 


2 


1168.90 


CIV 1168.993 


Ph 




1168.75 


CIV 1168.849 


Ph 


3 


1163.70 


7 


7 


4 


1158.97 


PII** 1159.087 


NN 


5 


1157.83 


CIUV15.01+16 


Nb,NN 




1157.72 


CIUV15.01+16 


Nb,NN 




1157.63 


CIUV15.01+16 


Nb,NN 




1157.48 


CIUV15.01+16 


Nb,NN 


6 


1155.73 


CIUV19 


Nb,NN 


7 


1152.71 


PII 1152.818 


NN 




1152.53 


PII 1152.818 


Nb 


8 


1148.11 


7 


7 


9 


1146.75 


O VI 1147.072 + 1146.791 


Ph 


10 


1144.85 


Fe II 1 144.938 


NN 


11 


1143.13 


Fe II 1 143.226 


NN 


12 


1142.28 


Fe II 1 142.366 


NN 


13 


1141.28 


7 


7 


14 


1139.50 


CI 1139.793 (UV22) 


Nb 


15 


1134.86 


NI 1134.980 


NN 




1134.70 


NI 1134.980 


Nb 




1134.30 


NI 1134.415 


NN 




1134.10 


NI 1134.415 + 1134.1653 


Nb+NN 




1133.90 


N I 1134.165 


Nb 


16 


1125.35 


Fell 1125.448 


NN 


17 


1124.65 


O VI 1124.809+1124.716 


Ph 


18 


1122.49 


O VI 1122.618 + CI UV27 


Ph,Nb 


19 


1122.25 


O VI 1122.348 


Ph,Nb 


20 


1121.88 


Fell 1121.975 


NN 


21 


1117.78 


C I UV29 


Nb.NN 




1117.44 


C I UV29 


Nb,NN 


22 


1107.48 


CIV 1107.593 


Ph 


23 


1106.21 


CI 1106.316 


NN 


24 


1096.78 


Fell 1096.877 


NN 


25 


1095.20 


7 


? 


26 


1085.42 


Nil** 1085.710 


Nb 


27 


1085.25 


Nil** 1085.551 + 1085.533 


Nb 


28 


1084.23 


Nil* 1084.584+ 1084.566 


Nb 


29 


1083.90 


Nil 1083.994 


NN 




1083.70 


Nil 1083.994 


Nb 


30 


1066.57 


Arl 1066.660 


NN 


31 


1063.09 


Fell 1063.176 


NN 


32 


1062.47 


S IV 1062.664 


Ph,Nb 




1062.35 


S IV 1062.664 


Nb 


33 


1048.14 


Arl 1048.220 


NN 




1047.98 


Arl 1048.220 


Nb 


34 


1039.12 


O I 1039.230 


NN 




1039.02 


O I 1039.230 


Nb 


35 


1037.48 


O VI 1037.617 


Ph 


36 


1036.75 


CII* 1037.018 


Nb 


36 


1036.08 


CII 1036.337 


Nb 


35 


1031.79 


O VI 1031.926 


Ph 


37 


1012.24 


S III 1012.495 


Nb 


38 


999.76 


7 


? 


39 


997.78 


P III 998.000 


Nb 


40 


991.38 


NIII 991.577 + 991.511 


Nb 




991.22 


NIII 991.577 + 991.511 


Nb? 


41 


989.79 


Si II 989.873 


NN 




989.56 


N III 989.799 


Nb 


42 


988.71 


O I 988.773 


NN 




988.47 


O I 988.773 + 988.6549 


Nb,NN 


43 


976.85 


C III 977.020 


Ph,Nb,NN 


44 


971.65 


O I 977.738 


NN 




971.48 


1 977.738 


Nb 


45 


964.54 


N I 964.626 


NN 




964.40 


N I 964.626 


Nb 




963.93 


N I 963.990 


NN 




963.74 


N I 963.990 +P II 963.801 


Nb 



TABLE 2 — Continued 



Feature 


A 


Atomic ID, ^„,,f 

A 


Component^ 


46 


954.02 


N I 954.104 


NN 




953.88 


N I 954.104 + 953.970 


Nb,NN 




953.73 


N I 953.970 


Nb 




953.60 


N I 953.655 


NN 




953.43 


N I 953.655 


Nb 




953.35 


N I 953.415 


NN 




953.21 


N I 953.415 


Nb 


47 


952.34 


N I 952.415 


NN 


48 


950.67 


1 950.885 


Nb 




950.47 


P IV 950.657 


Nb+(Ph?) 


49 


948.60 


O I 948.686 


NN 




948.47 


O I 948.686 


Nb 


50 


945.33 


C I 945.579 


Nb 




945.09 


C I 945.338 


Nb 




944.95 


CI 945.191 


Nb 


51 


944.40 


S VI 944.523 


Nb+Ph 


52 


940.01 


7 


7 


53 


936.55 


O I 936.623 


NN 




936.40 


O I 936.623 


Nb 


54 


934.85 


7 


? 




934.73 


7 


7 


55 


933.28 


S VI 933.378? 


7 


56 


931.90 


Aril* 932.054 


Nb 


57 


924.87 


1 924.950 


NN 




924.73 


1 924.950 


Nb 


58 


921.77 


1 921.857 


NN 




921.64 


O I 921.857 


Nb 


59 


920.11 


7 


7 


60 


919.85 


7 


7 


61 


919.54 


Ar II 919.781?,0 I 919.658? 


7 


62 


916.77 


O I 916.960+916.815 


Nb,NN 


63 


915.47 


Nil 915.61 


Nb,NN 



Ph = Photosphere, Nb = Nebular, NN = Non-Nebular, ? = Unknown 



The model atomic and molecular absorption spectrum was 
created by adding to a master optical depth template, cov- 
ering 900 - 1400 A, all the contributions from the atomic 
and molecular hydrogen lines for all the velocity components. 
There are two distinct components with heliocentric veloc- 
ities of -75 and -28 ± 2 km s^^ for molecular hydrogen. 
The first is associated with the nebular expansion and the 
second is associated with non-nebular ISM. For the atomic 
Lyman series we used transitions up to « = 49, with the ve- 
lo city offsets, column dens ities and doppler parameters listed 
in lMcCandliss et aP (l2007h . For hot mol ecular hydrogen we 
used the H200IS rovibrational templates of lMcCandlissI (l2003h 
with b = 7 km s^', offset bj^-75 km s^'and sca led to the 
column densities listed in iMcCandliss et alj ( l2007h . We ex- 
tended the calculation up to J"= 15 and v"= 3 by assuming 
the upper level ro-vibration populations were thermal with a 
temperature of 2040 K as given by the best fit single temper- 
ature model. We also used the H200IS templates to model the 
cold molecular component at -28 km s^' with a pure thermal 
distribution of 200 K. The transmission function, defined as 
the negative exponential of the master optical depth template, 
is multiplied to the stellar continuum model, which has been 
blueshifted to the systemic radial velocity of the nebula, -A2 
km s^\ to produce the model of continuum absorption by 
hydrogen. This model is also made available on the H200IS 
website. 

Figure |2a| shows the resulting continuum absorption model. 
The top spectral sequence in each panel is an overplot of the 
sl2 (orange) and s21 (green) spectral extractions with the ab- 
sorbed continuum model (black). The lower sequence in each 
panel is the residual formed by substracting the model from 
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Fig. 4.— Stellar photospheric O Vie and O Ve A1371.292 lines. Stellai- photospheric O Vie and O Ve >.1371.292 lines. The O Ve ^1371.292 line was used 
to established the gravitational redshift of the system. See Figure[3]for a description of the colors. 
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sl2 (orange) and s21 (green). An offset of 10 ergs cm 
s^' A^' and a two pixel smoothing has been applied to the 
residuals. The solid horizonal black line marks the zero level 
for the residuals. The high frequency spikes that go rapidly 
from positive to negative in the residual are created by lo- 
cal wavelength misalignment of the spectral extractions with 
respect to the molecular model. The strongest stellar photo- 
spheric, nebular and non-nebular absorptions by species other 
than H2, H I or photospheric He II are clearly revealed in the 
residuals. Emission lines of He II A,1084, A,992 and A,958 are 
also evident. The absorption systems are numbered in the fig- 
ure in reference to an entry in Table |2] where the line speci- 
fications are detailed. The identifications are not exhaustive; 
small unidentified lines may yet lurk in the noise. 

The residual spectrum allows us to assess the success of the 
absorbed continuum model in reproducing the features in the 
spectral extractions. Strictly speaking, the subtraction of fea- 
tures should take place in the optical depth space t(A,) as op- 
posed to the transmission space exp(— t(A,)). Abundance de- 
terminations whether by profile fitting or equivalent width de- 
termination must account for the optical depth blending cor- 
rectly. However, as our purpose here is merely line identifica- 
tion, the simple subtraction process is adequate. 

4. STIS AND FUSE ABSORPTION PROFILES AS A FUNCTION OF 
VELOCITY 

We separate the display of individual absorption profiles in 
Figures [3] -[TS] into four categories, stellar, nebular plus stel- 
lar, nebular CNO and nebular metals. The wider stellar lines 
we plot over a velocity bandpass of -250 to H-150 km s^'or 
-200 to H-lOO km s^^ Otherwise we use a bandpass of -150 
to H-50 km s^ ' . We omit the broad H I and He II photospheric 
features. The black dashed vertical lines mark the location 
of the hot molecular hydrogen component at -75 km s^' the 
systemic velocity at -42 km s^', and the cold molecular hy- 
drogen component at -30 km s^'. The red dashed vertical 
lines mark the gravitational redshift of the systemic velocity. 

The sl2 and s21 spectra are shown in orange and green 
respectively. The continuum model, including the hydro- 
gen lines, is plotted as a thin black line. Regions where the 
agreement is poor between the wavelength scales for the si 2, 
s21, and model spectra become readily apparent as do regions 
where the continuum flux is too high, especially near the O VI 
lines. For lines in the STIS region the spectra are plotted in 
green. 

4. 1 . Stellar Photospheric Features and the Absolute 
Wavelength Scale 

Comparison of the FUSE and STIS wavelength scales re- 
vealed a systematic offset. The reconciliation of this offset is 
essential for investigating the kinematics of the nebular out- 
flow, where we seek to determine the velocity of the various 
molecular and atomic features with respect to the systemic 
velocity of the nebula. In principle, lines that arise from the 
photosphere should match the systemic velocity of the system 
less the gravitational redshift offset. In practice this has been 
difficult to realize except for the narrowest of photospheric 
features. 

We have examined the overlap of the spectra in the com- 
mon wavelength regions below 1 190 A, where high excitation 
photospheric C IV and O VI lines are located. We have also 
checked the registration of the excited molecular hydrogen 
lines that appear in the STIS bandpass above 1 190 A with the 
continuum plus hydrogen model. The disagreements revealed 



from these examinations have been reconciled as described 
below. 

An absolute reference to the heliocentric velocity was es- 
tablished by close examination of the O V A,l 37 1.296 feature. 
This narrow line is a transition between two highly excited 
states in O V and is expected to be an excellent indicator of 
the photospheric restframe (Pierre Chayer private communi- 
cation) . For a compact object of the mass and radius found 
by Na piwotzkH (|l999) we expect the photospheric lines to ex- 
perience a gravitation redshift, Vgr — c{{l - ^^^yi — 1) = 

5.1 km s^'. Applying a shift of -13 km s^^ to the STIS spec- 
trum placed the centroid of the O V A ,l 37 1.296 at - 37 km s^ ' , 
as expected for a V^ys of -42 km s"' (IWilsonlll953h . 

In the original analysis of the FUSE M27 spectra by 
McCandliss (2001) the hot nebular molecular hydrogen com- 
ponent was defined to be at -69 km s^^ Examining the over- 
lap of common wavelength features in the FUSE and STIS 
spectra after defining the O V A-1 37 1.296 to be at -37 km s^\ 
we found it necessary to shift the FUSE spectra blueward 
by -6 km s^\ such that the hot nebular molecular hydro- 
gen is now at -75 km s^'. The most useful overlap lines 
for assessing the alignment were the doublet blend of C IV 
Ul 168.849, 1168.993 and the the narrow O VI Ul 171.56, 
1172.44. We note that the wavelengths of the O VI doublet 
given in the NIST online tables appear to be in error by « - 
0.42 A. lJahn et al.l(l2006h have produced an new empirical set 
of O VI wavelengths, which agree to within 0.04 A of those 
found here. We also examined the FUSE N I multiplets at 
Ul 134 - 1 135 and the STIS N I multiplets at U1200 - 1201 
to confirm the wavelength reconciliation. 

O Viand OVe— In Figure[3]and |4]we show the O VI and O V 
profiles we have identified as being photospheric in origin. 
Those lines that arise from absorption out of energy levels 
well above the ground state are designated as either O Vie or 
O Ve. The O VI A,1037.62 resonance line shows some very 
slight signs of blue shifted nebular absorption. 

ClVe — In Figure|5]we show the excited C IVe lines. None of 
these lines is saturated. Overlapping STIS spectra if available 
are plotted in purple (and grey if two orders are available). 
These spectra are comparatively noisy, but the alignment with 
the FUSE spectra agrees as well as the alignment between the 
sl2 and s21 spectra. We note that the error in the systemic 
velocity is of order the STIS resoution element and is three 
times the FUSE resolution element. We consider the agree- 
ment of line profiles from spectra acquired with two different 
instruments to be excellent. The STIS order, shown in grey, 
of the C IVe A.A.1 168.849, 1168.993 doublet has a spurious 
absorption feature that does not appear in the order shown in 
purple nor in the FUSE spectra and can be ignored. 

4.2. Photospheric + Nebular Features 

C IV and N V — The high ionization resonance doublets C IV 
U1548.204, 1550.781, and N V U1238.821, 1242.804 show 
signs of nebular absorption to the blue of -37 km s^ ' and pho- 
tospheric absorption to the red as can be seen in Figure|6l The 
nebular absorption component in the N V lines is strong only 
between -42 and -75 km s^' and is just barely saturated at 
« -60 km s^K In contrast, the nebular absorption component 
in the C IV lines spans -42 to -1 15 km s^' and is completely 
saturated from -50 to -95 km s^^ 
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Fig. 6.— Absoiption as a function of velocity for CIV XX1548.204, 1550.781, and N V XX1238.821, 1242.804, showing signs of nebular absorption blueward 
of -37 km s^' and photospheric absorption redward. See Figure[3]for a description of the colors. 



4.3. Nebular CNO 

In general, these lines show absorption from the nebula 
to the blue, and varying strengths of the intervening (non- 
nebular) ISM absorption components to the red of the sys- 
temic velocity. The nebular features that show absorption 
only in the nebular expansion (i.e. at velocities blueward of 
the systemic velocity) are from intermediate or low ionization 
species. The non-nebular velocity features are most promi- 
nent in the lowest ionization and neutral species along the 
line-of-sight. 

cm and C//— We show Hnes of C III and C II in Figure |2l 
Like the C IV lines, these lines are heavily saturated through- 
out the nebular flow region blueward of —42 km s^^ The sat- 
uration makes it difficult to tell, which if any, ion is dominant 
in the different nebular velocity regimes. 

N III, N II and N 1 — The N III lines in Figure [8] are strongly 
blended with overlapping molecular hydrogen. The blended 
profile indicates the absorption is saturated between -60 and 
-100 km s^'. The N II* and N II** lines in Figure|9]are sim- 
ilarily messy. However, N II A,1083.994 is relatively clean. 
It shows absorption throughout the flow, being less saturated 
at low velocities and becoming completely saturated at -75 



km s^'. This line also shows stonger nebular absorption than 
non-nebular absorption emphasizing its dominance over N I, 
which shows the opposite behavior The N I multiplets at 
U1134.165 - 1134.980 and Ul 199.550- 1200.710, in Fig- 
ure [TOal simply show unsaturated absorption centered at -75 
km s^'. 

01 — The O I lines in Figure [TT] show a range of saturated 
to unsaturated profiles centered at -75 km s^'. It may be 
possible to determine the column density of this species fairly 
accurately with a curve of growth, after properly accounting 
for the continuum placement and molecular hydrogen optical 
depth subtraction. 

4.3.1. Nebular Metals 

The transition zone between high ionization and low ion- 
ization occurs at the velocity of -75 km s^^, where H I, C I, 
N I, O I and molecular hydrogen show up most strongly in the 
nebula. 

Si IV, Si III and Si II — The high ionization low velocity, low 
ionization high velocity dichotomy is most clearly seen in the 
Si IV X1393.760 and Si II A.1 193.290 lines as shown in Fig- 
ure [T2a] Since these lines have nearly equal log(/A-) they are 
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Fig. 7. — Absoiption as a function of velocity for the C III and C II lines, 
region blueward of —42 km s" ' . See Figure[3]for a description of the colors. 



Like the C IV lines, these lines ai'e heavily saturated throughout the nebular flow 
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lines are also blended with each other. See Figure |3]for a description of the 
colors. 



reliable indicators of the relative ionization as a function of 
velocity. Si II is stronger than Si IV in the zone between -75 
and -1 10 km s^', while Si IV is stronger than Si II in the zone 
between -A2 and -75 km s^^. The Si III A.1206.500 line is 
saturated throughout most of the flow. It is a very strong line 
with a log(/A-) = 3.29, which makes it difficult to determine if 
this species is the dominant ion throughout the flow. 

P IV, P III and P II— We detect weak P IV A.950.657, P III 
A.998.000, and P II A.1 152.818 at progressively high velocities 
of « -65, -70, -80 km s^' respectively as shown in Figure[T3] 
The P V A.1 1 17.977 line (not shown) of the P V doublet is very 
weak, while A.1 128.008 is blended with H2. 

S IV, S 111 and 5 // — The S IV - S II ions show similar behavior 
to Si and P in Figure [H The S III Al 190.203 fine appears to 
have a feature near -42 km s^'. However, this is just the 
nebular component of the nearby Si II Al 190.416 line. The 
nebular portion of this feature is clean and agrees well with 
the shape of the S III A1012.495 profile, with the absorption 
stongest near -70 km s^'and weakening slowly to the blue 
and relatively quickly to the red. 



Ar II and Ar I— Figure [15] shows Ar II A919.781, Ar II* 
932.054, and Ar I 1048.220 are detected, with Ar II appear- 
ing at w -55 km s^' and Ar I appearing at « -75 km s^^ 
Ar I 1066.660 is blended with an overlapping H2 line. Ar II 
919.781 is blended with a nearby molecular hydrogen line but 
the nebular Ar II* A932.054 is clearly visible without model 
subtraction. The thermal production of Ar II* requires a tem- 
perature ^ 10,000 K, and is a confirmation that the gas at low 
velocities is much hotter than that indicated by the molecular 
hydrogen. The absorption of the Ar II lines is stronger than 
the Ar I lines even though they have transition strengths a fac- 
tor of 10 lower This suggests the Ar II is the dominant over 
Ar lin the nebula. 

Fe III and Fe II~ Fe III Al 122.524 may be blended with a 
molecular hydrogen line aon the blue side of the photospheric 
O Vie AAl 122.593 - 1 122.618 doublet as shown in Figure[T6] 
However, without a good idea of what the photospheric line 
shape is in this region, we cannot claim a detection. The 
strongest transitions of Fe II A1608.451 and Al 144.938 appear 
weakly at -75 km s^'. These lines have transition strengths 
(log (/A) = 2.06, 2.11) that are nearly identical to that of Si II 
A1304.370 (log (/A) = 2.08), yet the Si II feature is much 
stronger The ionization potential of Si II (16.3 eV) is also 
nearly identical to that of Fe II (16.2 eV) a nd th ey also have 
similar condensation temperatures (iMortonll2003() . We con- 
clude that Fe II relative to Si II is underabundant in the nebula. 

Arguing in a similar vein, we see that the transition at 
Fe II Al 133.665 is not detected in the nebula at all. It has 
nearly the same transition strength (log (/A) = 1.15) as the 
S II A1253.81 1 line (log (/A) = 1.14), which is weak but eas- 
ily seen. We conclude that Fe II relative to S II is also un- 
derabundant. The solar abundance of Si and Fe are nearly 
identical (7.56, 7.50 respectively) and are slightly higher than 
S (7.2). It is not unusual to have the Si abundance greater 
than the Fe in the typ ical ISM because it is depleted onto dust 
jSavage & Sembadi|[T996i) . However, here we have no Fe in 
the gas phase and appearently very little dust in the diffuse 
nebular medium. This is a puzzle. 



5 . SUGGESTIONS FOR FUTURE INVESTIGATIONS 
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Fig. 9. — Absorption as a function of velocity for the N II >.^1083.994 - 1085.710 multiplet. The N II* and N II** lines ai'e blended with each other and 
molecular hydrogen. The N II line is relatively clean. See Figure[3]for a description of the colors. 
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Fig. 10a. — Absorption as a function of velocity for the N I lines show unsaturated absorption at -75 km s ' . See Figure[3]for a description of the colors. 




Fig . 1 Ob . — More N I lines. See Figure[3]for a description of the colors. 



M27 is a excellent laboratory for testing theories regarding 
the abundance kinematics in PNe because the stellar temper- 
ature, mass, gravity, and the nebular mass, distance, abun- 
dances, extinction and excitation states, are well quantified. 
The reduced spectra provided by this study and the associated 
atomic and molecular hydrogen model for the nebular and 
non-nebular absorptions, should enable a photospheric analy- 
sis of the metal abun dances, similar to the effort undertaken 
by lJahnet all (l2006l) for PG 1159 - 035. Establishing a re- 
liable stellar continuum is the first step towards determining 
accurate nebular absorption line abundances for comparision 
with those derived from the emission line analyses of lBarkerl 
([1984) and Hawley & Miller ( 1978). 

The finding in Paper I, of little extinction by dust in the neb- 
ula raises a number of interesting questions concerning the 
depletion of metals in the diffuse and clumpy media, which 
go beyond the scope of this investigation. A major limitation 
to this effort is the requirement for a good stellar model that 
can reliably reproduce the observed photospheric features and 
continuum, especially in the vicinity of the blend of nebular 
Fe III A.1 122.52 with photospheric O VI A.1 122.62. If such a 
stellar model can be produced, the total S, Si, P and Fe abun- 
dances should be derivable by profile modeling after account- 
ing for the atomic and molecular hydrogen absorption pro- 
vided here. The question of whether the abundances of these 
metals change across the neutral transition zone could pro- 
vide information on whether photo-evaporation of the glob- 
ules is an important process for the enrichment of metals in 



the high velocity zone. A model of the nebular ionization as a 
function of velocity would be useful for this purpose. Such a 
model might be produced by melding the detailed time depen- 
dent approach to the ca lculations of the atomic and molecular 
emissions followed by iNatta & HoUenbachl (I 1998'). with the 
radiative hydrodynamical rigor of FVillaver et a l. (2002). The 
problem may also be approached by us ing the Sobolev plus 
exact integration (SEI) method of Lamers et al.l ([T987) . from 
which a empirical parameterization of the wind velocity 'law" 
could be obtained. 
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Fig. 12a. — Compai'ision of the Si IV, Si III and Si II. Si II dominates Si IV in the zone between -75 and -110 km s ' while Si IV dominates the lower 
ionization species near -60 km s^ ' . See Figure[3]for a description of the colors. 
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Fig. 13.— P IV W50.657, P III A998.000, and P II U 152.818 lines aie 
weak. Although the P IV and P III are surrounded by molecular hydrogen 
they show behavior similar to Si. See Figure|3]for a description of the colors. 
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Fig. 14. — S IV - S II lines show behavior similai' to Si. See Figure|3]for a description of the colors. 
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Fig. 16. — The detection of Fe III cannot be claimed. Fe II is underabundant relative to S II. See Figure[3]for a description of the colors. 



